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Abstract 

We prove that isospectral simple analytic surfaces of revolution are isometric. 

Introduction 

This article is concerned with the inverse spectral problem for metrics of revolution on S 2 . We will assume 
that our metrics are real analytic and belong to a class TZ* of rotationally invariant metrics which are of 
'simple type' and which satisfy some generic non-degeneracy conditions (see Definition (0.1)). In particular, 
we will assume they satisfy the generalized 'simple length spectrum' condition that the length functional on 
the loop space is a clean Bott-Morse function which takes on distinct values on distinct components of its 
critical set (up to orientation). Denoting by Spec(S 2 7 g) the spectrum of the Laplacian A g , our main result 
is the following: 

Theorem I Spec: TZ* -> R +IN is 1-1. 

Thus, if (S 2 ,g), (S 2 ,h) are isospectral surfaces of revolution in TZ* , then g is isometric to h. It would 
be very desirable to strengthen this result by removing the assumption that h e TZ* , thereby showing 
that metrics in TZ* are spectrally determined within the entire class of analytic metrics on S 2 with simple 
length spectra. The only metric on S 2 presently known to be spectrally determined in this sense is the 
standard one (which is known to be spectrally determined among all C°° metrics). A metric h satisfying 
Spec(S 2 , h) — Spec(S 2 , g) for some g e TZ* must have many properties in common with a surface of revolution 
of simple type; it would be interesting to explore whether it must necessarily be one. 

Let us now be more precise about the hypotheses. First, we will assume that there is an effective action 
of S 1 by isometrics of (S 2 ,g). The two fixed points will be denoted N, S and (r, 9) will denote geodesic polar 
coordinates centered at N , with 9 = some fixed meridian 7m from N to S. The metric may then be written 
in the form g = dr 2 + a{r) 2 d9 2 where a : [0, L] -> IR+ is defined by a(r) = ±\S r (N)\, with \S r (N)\ the 
length of the distance circle of radius r centered at N. For any smooth surface of revolution, the function a 
satisfies a (2p) (0) = a (2p) (L) = 0,a'(0) = l,a'(L) = -1 and two such surfaces (S 2 ,gi) (i = 1,2) are isometric 
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if and only if L\ = L 2 and ai(r) = a 2 (r) or ai(r) = a 2 (L — r). We will then assume that the metrics belong 
to the following class 1Z of simple analytic surfaces of revolution [Bl] : 

(0.1) Definition 1Z is the moduli space of metrics of revolution (S 2 , g) with the properties: 

(i) g (equivalently a) is real analytic; 

(ii) a has precisely one non-degenerate critical point r <E (0, L), with a"(r ) < 0, corresponding to an 
'equatorial geodesic' je; 

(Hi) the (non-linear) Poincare map V JE for ^e is of twist type (cf. §1). 

We denote by 1Z* C 1Z the subset of metrics with 'simple length spectra' in the sense above. 

Regarding 'simple length spectra,' we recall that the closed geodesies of a surface of revolution come in 
one-parameter families of a common length, filling out invariant torii T for the geodesic flow. The canonical 
involution a(x, £) = (x,—£) takes T to its 'time reversal' — T, and takes the closed geodesies of T to 
their reversals on — T. A closed geodesic and its reversal have the same length, so the length spectrum is 
automatically double except for the length 2L of the torus of meridians, which is cr-invariant. The simple 
length spectrum hypothesis is that up to time reversal, the common lengths of the closed geodesies on distinct 
torii are distinct (cf. Definition 1.2.2). In fact, it would be sufficient for the proof that the length 2L of the 
'meridian' closed geodesies is not the length of closed geodesies on any other torus. In any case, it is not 
hard to show that 7Z* is residual in 1Z (cf. Proposition 1.2.4). 

The condition (iii) actually appears in the proof in the following way: the quadratic coefficient 'a := h"(0)' 
of the classical Birkhoff normal form of the metric |£| g at the torus of meridian geodesies must be non- 
vanishing (see Definition (1.4.5) for the precise meaning of /i(£))- This condition is used in Proposition 
(4.1.2) and Corollary (4.1.3) to evaluate the wave invariants for the meridian torus. 

As will be explained further below, the main purpose of the non-degeneracy and simplicity conditions 
is to insure that there are global action-angle variables for the geodesic flow. These conditions rule out 
several types of surfaces of revolution: First, they rule out Zoll surfaces of revolution, which are degenerate 
in every possible sense. It is indeed unknown at this time whether real analytic Zoll surfaces of revolution 
are determined by their spectra. They also rule out 'peanuts of revolution' (which have hyperbolic waists) 
and other natural rotational surfaces such as Liouville torii. 

To our knowledge, the strongest prior result on the inverse spectral problem for surfaces of revolution is 
that of Bruning-Hcintzc [B.H]: smooth surfaces of revolution with a mirror symmetry thru the x — y plane 
are spectrally determined among metrics of this kind. There are also a number of proofs that a surface of 
revolution is determined by the joint spectrum of A and of d/d6, the generator of the rotational symmetry 
[Kac][B][Gur]. 

The method-of-proof of Bruning-Hcintzc was based on the observation that the invariant spectrum can 
be heard from the entire spectrum. Hence by separating variables the problem can be reduced to the inverse 
spectral problem for ID even singular Sturm-Liouville operators, which was solved by Gelfand-Levitan and 
Marchcnko. 

Our proof of the Main Theorem is based on different kind of method, and the inverse result presented 
here is hopefully just one illustration of it. We begin with the facts that simple surfaces of revolution are 
completely integrable on both the classical and quantum levels and that the Laplacian has a global quantum 
normal form in terms of action operators. We then study the trace of the wave group and prove that from 







its singularity expansion we can reconstruct the global quantum normal form. Finally we show that this 
normal form determines the metric. 

This approach is suggested by the recent inverse result of Guillemin [G.l], which shows that the microfocal 
normal form of A around each non-degenerate elliptic closed geodesic can be determined from the wave 
invariants (see also [Z.1,2]). However, there is but one non-degenerate closed geodesic on a simple surface 
of revolution, so the direct application of this inverse result does not take full advantage of the situation. 
Rather, it is natural to start from the fact that the wave group is completely integrable in the following 
strong sense: namely, it commutes with an effective action of the torus S 1 x S 1 by Fourier Integral operators 
on L 2 (S 2 ). That is, there exist global action operators I\, I2 and a polyhomogeneous symbol H of degree 1 
on 1R 2 — such that y/~K = H(Ii,l2)- This is the global quantum Birkhoff normal form alluded to above. 
Our principal tool is the following inverse result: 

Main Lemma The wave trace invariants of (S 2 ,g) with g £ TV determine the quantum normal form H. 

This Lemma does not actually require that (S 2 , g) be a surface of revolution, but only that the geodesic 
flow is toric integrable, i.e. commutes with an effective Hamiltonian torus action. It immediately implies 
that the principal symbol H{I\, I2) of H(I\,l2) is a spectral invariant. Since H(Ii,I 2 ) is essentially a global 
Birkhoff normal form for the metric, the wave invariants determine the symplectic equivalence of the geodesic 
flow. Thus we have: 

Corollary 1 From the wave trace invariants of (S 2 ,g) with g £ TV we can determine the symplectic 
equivalence class of G l . 

Corollary 1 does not however determine the isometry class of a general g £ TV: As will be discussed in §1 
(see also [C.K]), simple surfaces of revolution are not symplectically rigid unless they arc mirror symmetric. 
Otherwise put, recovery of the classical Birkhoff normal form only determines the even part of the metric in 
the following sense: 

Corollary 2 A metric g £ TZ may be written in the form g = [f (cosu)] 2 du 2 + [sinu] 2 d9 2 (§1, [Besse]). 
From the wave invariants one can determine the even part of f. In particular if g is mirror symmetric, one 
can determine g among mirror symmetric metrics in TZ* from its spectrum. 

It is interesting to note that symplectic rigidity in the mirror-symmetric case gives a new and self- 
contained proof of the Bruning-Heintze result, without the use of Marchcnko's inverse spectral theorem for 
singular Sturm Liouville operators. Although this result is superceded by the Theorem, it may have some 
future relevance to other inverse problems. 

To complete the proof of the Theorem, we therefore have to study the subprincipal terms in H. The 
result is: 

Final Lemma From H one can determine the isometry class of g. 

It is in this last Lemma that we use in full that (S 2 , g) is a surface of revolution rather than just a surface 
with toric integrable geodesic flow. We also use in full that H is a global quantum normal form rather than 
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a microlocal one at je- In subsequent work we will investigate the analogues of the Final Lemma for the 
microlocal normal form at 7b, for more general toric intcgrable metrics and for metrics isospectral to toric 
integrable ones. 

To close this introduction, we discuss some background and some open problems connected with this 
work. 

The principal motivation for studying the inverse spectral problem for surfaces of revolution is its sim- 
plicity. There are to date very few inverse results which determine a metric from its Laplace spectrum within 
the entire class of metrics, or even within concrete infinite dimensional families. To our knowledge, only the 
standard S n for n < 6 and flat 2-torii are known to be spectrally determined. Hence it is desirable to have 
a simple model of how an inverse result might go. 

A second motivation is a somewhat loose analogy between surfaces of revolution and planar domains. 
Namely, in both cases the unknown is a function of one variable (the profile curve, resp. the boundary) 
which completely determines the first return times and angles of geodesies emanating from a transversal. 
That is, paths of billiard trajectories (broken geodesies) on a bounded planar domain are determined from 
collisions with the boundary, while paths of geodesies on surfaces of revolution are determined from collisions 
with a meridian (or with the equator). In analytic cases, it is plausible that the unknown function may be 
determined in large part by the spectrum of first return times from the local transversal. In the case of an 
analytic surface of revolution, this length spectrum determines the corresonding Birkhoff normal form and 
hence by Corollary 2 it determines the even part of the profile curve. Similarly, in the case of an analytic 
plane domain, it is proved in [CV.2] that the even part of the boundary of an analytic domain may be 
determined from the Birkhoff normal form of the billiard map at a bouncing ball orbit. Hence, there is a 
similarity in the relation between the unknown function and the local classical Birkhoff normal forms. It is 
interesting to observe in this context that the rigidity result of Colin de Verdiere [CV.2] is quite analogous 
to the inverse result of Bruning-Hcintze. 

Some immediate open problems: First, there is the symplectic conjugacy problem mentioned above. 
Second, can one relax analyticity to smoothness in the Theorem above? This is likely to follow from a more 
intensive analysis of the wave invariants. Third, can one extend it to other 'non-simple' types of surfaces 
of revolution? The main obstacle is that one will generally not have global action-angle variables or global 
quantum normal forms. What about completely integrable systems in higher dimensions? 

Finally, we would like to thank D.Kosygin for several helpful conversations on this paper, and B.Kleiner 
for giving us up to date information about the status of the conjugacy problem for geodesic flows on surfaces 
of revolution. 

1 Classical dynamics 

1.1 Global action-angle variables 

From a geometric (or dynamical) point of view, the principal virtue of metrics in 1Z is described by the 
following: 

(1.1.1) Proposition Suppose g is a real analytic metric of revolution on S 2 such that a has precisely 
one non-degenerate critical point at some r € (0, L). Then the Hamiltonian \(\ g := y/Yl 9 i: *€i€j on T*S 2 is 
completely integrable and possesses global real analytic action-angle variables. 
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Proof: The complete integrability of |£| g (i.e. of the geodesic flow) is classical, and follows from the 
existence of the Clairaut integral pe(v) := (v, Since the Poisson bracket {pg, \£\ g } = 0, the geodesies 
are constrained to lie on the level sets of pg; and since both |£| g and pg are homogeneous of degree one, the 
behaviour of the geodesic flow is determined by its restriction to S*S 2 = = 1}. With the assumption 
on a, the level sets are compact and the only critical level is that of the equatorial geodesies 7^ C S*S 2 
(traversed with either orientation). The other level sets are well-known to consist of two-dimensional torii. 
(Had we allowed the existence of at least two critical points in a, there would exist a saddle level, i.e. an 
embedded non-compact cylinder). 

The existence of global action-angle variables follows from the general results of [D][G.S] and have been 
constructed explicitly for simple surfaces of revolution in [CV.l]. The general formula is as follows: Let 

(1.1.2) P = m g , P g) : T*S 2 —>£?:= {(61,62) : N < «(r )6i} C M x IR+ 

be the moment map of the Hamiltonian H 2 -action defined by the geodesic flow and by rotation. The singular 
set of P is the closed conic set Z := {(r a ,9, 0,p$) : 6 e [0, 2n), pg e M}, i.e. Z is the cone thru the equatorial 
geodesic (in either orientation). The image of Z is the boundary of B; the map P\T'S g s 2 -z 1S a trivial 
S 1 x S 1 bundle over the open convex cone B Q (the interior of B). For each b e B Q , let Hi(Fi,,ZZ) denote 
the homology of the fiber Fb := P _1 (6). This lattice bundle is trivial since B is contractible, so there exists 
a smoothly varying homology basis {71(6), 72(6)} € Hi(Fb, ZZ) which equals the unit cocircle S* N S 2 together 
with the fixed closed meridian 7m when b is on the center line 1R + • (1, 0). The action variables are given by 
[CV.l, §6] 

r f 1 rr+(b) I l2 

(1.1.3) h(b)= Zdx=pg, h{b) = £dx=- bl-^dr + \b 2 \ 

where r±(b) are the extremal values of r on the annulus 7r(Ffc) (with it : S*S 2 — > S 2 the standard projection). 
On the torus of meridians in S*S 2 , the value of I2 equals ^ and it equals one on the equatorial geodesic. So 
extended, I\, I2 are smooth homogeneous functions of degree 1 on T*5 2 , and generate 27r-periodic Hamilton 
flows. 

It follows that the pair 1 := (hjh) generates a global Hamiltonian torus (S 1 x S' 1 )-action commuting 
with the geodesic flow. The singular set of 1 equals Z := {I 2 = ±pg}, corresponding to the equatorial 
geodesies. The map 

(1.1.4) I : T*S 2 -Z^T := {(x, y) e JR. x IR+ : \x\ < y} 

is a trivial torus fibration. Henceforth we often write Tj for the torus I _1 (7) with I e r o and let T denote the 
closure of r o as a convex cone. The symplectically dual angle variables (<pi = 9, (j) 2 ) then give, by definition, 
the flow times (mod 2tt) along the orbits of (Ii,h) from a fixed point on which we may take to be the 
unique point lying above the intersection of the equator and the fixed meridian on Ff, with the geodesic 
pointing into the northern hemisphere. 

So far, we have only assumed the metric to be C°°. We now observe that if g is real analytic, then so are 
Ii, J2. This is obvious in the case of I\ and follows from the explicit formula (1.1.3) for 12- ■ 

Since the metric norm function |£| g commutes with I\,I 2 , it may be expressed as a function H(I) of the 
action variables. Hamilton's equations for the geodesic flow then take the form 

4 = 0, 0fc=w fc (J), (fc = l,2) 
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where 

(1.1.5) w(/)=V/i?(7) 

is the frequency vector of the torus Tj with action coordinates /. The geodesic flow on Tj is then given by 

(1.1.5a) G t (I,4>) = (I,(f> + tu;i) 

so that all the geodesies are quasi-periodic in action-angle coordinates. 

The frequency vector u>j is homogeneous of degree on T*5 2 - 0, and hence is constant on rays of torii 
IR + T/. To break the IR + symmetry we restrict to the level set {H(I) = 1} C T in action space and view 
the frequency vector as the map: 

w : {H = 1}^R 2 . 

Since ViH(I) _L Tj({H = 1}) the frequency map is more or less the Gauss map of {H = 1} (although 
it is not normalized to be of unit length). As a map of the global action space, the frequency map is the 
Legendre transform associated to H (cf. [F.G, p. 338]). 

(1.1.6) Definition We say that the simple surface of revolution (S 2 ,g) is globally non-degenerate if 
u)\{h=i} is an embedding. 

This is the natural homogeneous analogue of the non-degeneracy condition of [F.G, loc.cit] to the effect 
that the Legendre transform be a global diffeomorphism, and has previously been studied in some detail by 
Bleher in the setting of simple surfaces of revolution [Bl,§6]. As will be seen in §1.3, the curve {H — 1} is 
the graph of a smooth function of the form I2 — F(Ii) in the cone r o , and the non-degeneracy condition 
(1.1.6) will follow as long as F is a convex or concave function. In fact, in the proof of the Theorem we will 
only need to use that {H = 1} is non-degenerate at the one point I2) — (0, 1). This is sufficient because 
we assume the metric to be real analytic. 



1.2 Length spectrum and periodic torii 

We now come to the definition of length spectrum and simple length spectrum for a completely integrable 
geodesic flow. We first observe that the orbit thru (I, <j)) is periodic of period L if and only if 



(1.2.1a) Lwj = M € ZZ 2 

for some M 7^ 0. The minimal positive such L will be called the primitive period; the corresponding M is 
known as the vector of winding numbers of the torus Tj. M parametrizes the homology class of the closed 
orbit 7 since the latter has the form Y^j=i relative to the homology basis 

Due to the homogeneity, the period and vector of winding numbers are constant on the ray JRTj. By 
Eulcr's formula we then have 

ViH ■I = w I -I = H 
hence the length is given in terms of the winding vector by 
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Figure 1: 



or simply L = M ■ I on the unit tangent bundle H — 1. 

It is clear that the periodicity condition Luj — M is independent of cf>. Hence, all of the geodesies on Tj 
are closed if any of them are. ( This also follows, of course, from the transitivity of the torus action on each 
invariant torus.) We therefore say: 

(1.2.2) Definition 

(a) A torus Tj is a periodic torus if all the geodesies on it are closed. 

(b ) The period L of the periodic torus is then the common period of its closed geodesies. 

(c) The length spectrum C of the completely integrable system is the set of these lengths. 

(d) The completely integrable system has a simple length spectrum if there exist a unique periodic torus (up 
to time reversal) of each length L G C. 

In the last statement (d) we are referring to the canonical involution a : (x,£,) — ► (x, — £), which reverses 
the orientation of the geodesies. It is obvious that if Tj is a periodic torus of period L, then so is <r(T/). 

Let Per C S*S 2 — denote the set of periodic points for the geodesic flow on S*S 2 , i.e. the set of points 
which lie on a closed geodesic. It is a union of periodic torii in S*S 2 together with points along the equatorial 
geodesies (which are degenerate torii). The set of all periodic points in T*S 2 — is then equal to IR + Per. 
Since the invariant torii are parametrized by the points I e T of action space, it is convenient to parametrize 
Per by a subset of the level set {H(I) = 1} c T. 

(1.2.3) Definition TTie set of points I £ {H = 1} C T such that Tj C Per will be called, with a slight abuse 
of notation, the set of periodic points on {H — 1} and will be denoted by V . That is, V = {I G {H = 1} : 
3L e JR + ,Llu! g 2Z 2 }. 
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The following proposition will be needed later on (Proposition (4.1.4). 

(1.2.4) Proposition If(S 2 ,g) is non-degenerate (1.1.6), thenV is dense in {H = 1}. 

Proof: Let Q := uj(V) equal the image of V under the frequency map w. Then by definition, Q (for 'rational 
points') is the projection to the curve lu({H = 1}) of the integer lattice in ZZ 2 . It is clear that Q is a dense 
set in uj({H = 1}); since w is an embedding, V is dense in {H = 1}. ■ 

The next proposition shows that our inverse result is valid for a residual set of simple analytic surfaces 
of revolution. 

(1.2.5) Proposition Let 1Z* C 1Z be the subset of metrics with simple length spectra. Then 1Z* is a residual 
subset oflZ. 

Proof: If L, V E C, then there exist M, M' e ZZ 2 and 1,1' € Q such that L = M-I,L' = M' ■ I 1 . So L = V 
implies that M ■ I — M' ■ I 1 = hence that the /-coordinates of /, /' are dependent over the rationals. Since 
the length spectrum moves continuously and non-trivially under deformations in 1Z, such a dependence for 
fixed M, M' can only hold on a closed nowhere dense set. The proposition follows. ■ 

1.3 First return times and angles 

Let us consider more carefully the geometric interpretation of the Clairaut integral on a torus Tj C S*S 2 . 
Since I\ and H = |£| are independent commuting coordinates, and since there are global action-angle 
variables, the different invariant torii in S*S 2 are parametrized by the values I\ = i € [— 1, 1] of the Clairaut 
integral along {H = 1}. Indeed by (1.1.3), the second action coordinate I2 is determined from Ji on {H = 1} 
by the formula 

1 f-r+ih) I r2 

(1.3.1) I 2 = F(h) := |/!| + - / Jl- -^dr 

where r±(Ii) are the roots of I 2 = a(r) 2 . The torus in S*S 2 with Ii = 1 is therefore T( ti F(i)), which we will 
denote simply by T L . The projection of T L to S 2 is an annulus of the form r+(t) < r < r_(i). the geodesies 
on T/ project to S 2 as almost periodic curves oscillating between the two extremal parallels r = r + (t) and 
r = r_(i). 

We observe then that T L contains a unique geodesic r y l which passes thru the intersection of the reference 
geodesies 7m and 7^; 1 equals a(r )cosa where a{i) is the angle between j t and 7m- In other words, a(t) is 
the common angle with which the geodesies on T L intersect the meridians as they pass thru the equator in 
the direction of the northern hemisphere. Since the length Le of the equator 7^ equals 27ry / a(r c ), and since 
this length is a symplectic invariant of the geodesic flow, the coordinates 1 and a are related in an essentially 
universal fashion. Hence, cither I\ or a e (0, n) could be used as an action coordinate on {H = 1}; a is 
perhaps more geometric but I\ is more convenient in calculations. 

The picture is the same for any invariant torus Tj: Under the H + action on T* S 2 — 0, it scales to 
a torus T L C S* S 2 and all features of its geometry are identical to that of T L . Thus it carries a unique 
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(parametrized) geodesic 7/ such that 7/(0) is at the intersection je n 7a/, and so that 7j(0) points to the 
northern hemisphere. The initial angle variables of (77(0), 7/(0)) are therefore (</>i(0), ^2(0)) = (0,0). At 
time t we denote the angle variables of (ji(t), 7/(i)) by (<^i(i), 4>2{t)) where </>i essentially measures the 
meridianal angle and </>2 measures the equatorial angle. We then introduce the following 'first return times': 

(1.3.2) Definition 

(Ei) The equatorial first return time is the minimal time te(I) > such that 4>i{te{I)) = 2ir; 

(Eii) The equatorial first return angle loe(I) '■= 4>2(te(I)) is the change in angle along the equator of a 

geodesic on Tj leaving 775 fl 7 a/ at t — 0, upon its first return time to 7/5; 

(Mi) The meridianal first return time is the minimal time tm(I) such that 4>2{tm(I)) = 27r. 

(Mii) The meridianal first return angle is the angle change lom{I) '■= 4>i{ t m{I)) along the meridian of a 

geodesic on 7a/ leaving je H 7a/ at t = 0, upon its first return time to 7a/. 

The terminology 'first return time' is taken from dynamics. Note that <f>i = 0(mod2ir) is the equation of 
the curve on Tj which lies over the equator. Hence, Te(I) is the time of first return of ji(t) to the equator 
in the direction of the northern hemisphere. This is actually the second time of intersection of 7/ with the 
equator, the first one occurring when 77 (T) is heading to the southern hemisphere. This intersection is not 
in the projection of <pi = Q{mod2ir). Similarly, 4>2 = 0(mod2Tt) is the equation of the fixed meridian 7a/, and 
so ta/ is the time of first return to the arc jm (half of the closed geodesic). 
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Figure 2: 



The following gives some relations between the various angles and return times. 



(1.3.3) Proposition Let ujj = (101,1^2) be the frequency vector of the invariant torus Tj. Then: 
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(a) t e u) 2 = uj e ; 

(b) T M OJl = lu m ; 

(°) 2lT Z^T) = w m, 27T^g = UJ E . 

Proof: The equation of the geodesic (7/(£), 7/(i)) on Tj beginning at (0i,02) in angle variables is (0i + 
tu>i{I), + tLo 2 {I))- When the torus is projected to the base, then 0! measures the meridianal angle and 2 
measures the equatorial angle. So by definition of return times, 

(01 + Wl, 02 + T E U) 2 ) = (01 + 2-7T, 2 + UJe) 
(01 +T M Ul,<h +TMU 2 ) = (01 + WM,02 +27T). 

The statements in the proposition follow immediately. ■ 

These first return times (and angles) arc closely related to the (non-linear) Poincare maps of the geodesic 
flow. We recall that for each closed geodesic 7, the Poincare map "P 7 is defined as the first return map of 
the geodesic flow, restricted to a symplectic transversal S-y C S*S 2 (surface of section). It is well-known 
is a symplectic map [K]. Since the symplectic form on a cotangent bundle equals da (with a 
the action form), it follows that V*{a) — a is a closed 1-form on S 1 . Since 5 7 may be assumed contractiblc, 
it follows that there exists a function r 7 such that [P*(a) — a]\s y = dr 7 . Since the integral of a over an arc 
of a unit speed geodesic just gives its length, r 7 is the first return time of geodesies near 7 to 5 7 . 

In particular, let 7 = -f E be the equator (in one of its orientations), and let S E denote a symplectic 
transversal at the point (7b(0), 7^(0)) thru the fixed meridian jm- Since jm is transverse to the equator, 
we may define S E to consist of a small variation of 7^(0) moved up and down a small arc of jm- We see 
then that the first return time r 7E is precisely the first return time tm defined above. We also observe that 
the foliation of S* S 2 by invariant torii restricts to a foliation of S E by invariant circles for V lE , closing in 
on the center point where j E intersects S E - As noted above, the action coordinate I\ — l gives a natural 
action (radial) coordinate on S* 7 . Since 5 7 C S*S 2 , the H + homogeneity is broken and we may reformulate 
the twist condition (0.1 (iii)) as follows: 

(1.3.4) Proposition The Poincare map V 1E is a twist map of S E iff to' M 7^ 0, where oj' m — §i^jm- 

Proof: The coordinates (7i = i, <f>i) restrict to a system of symplectic coordinates on S E} in terms of which 
the Poincare map takes the form 

V 1B {i,(j)l) = (L, 01 +T M UJl) = (L, 01 + Wtf). 

By definition it is a twist map if oj' m 7^ in a neighborhood of t = in S E . m 
For background on the twist condition in a related context see [F.G][P]. 

The situation for the other periodic orbits is different since they come in one-parameter families. Thus, 
for the (closed) meridian geodesic 7m , a transversal Sm is given by the equator j E and a small variation of 
7^(0) along it. The foliation by invariant torii restricts to Sm to a foliation by invariant lines (non-closed 
curves), including the curve of closed geodesies thru -fM- The Poincare map P lM is then of parabolic type; 
the first return time is t e above. 
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For the inverse problem it will be necessary to have expressions for these return times and angles in 
terms of the metric. This will also make the twist condition more transparent. For ease of quotation, it is 
convenient to make a change of dependent and independent variables, following [Besse] and Darboux [D]. 
Equivalent expressions in the original polar coordinates can also be easily derived, and will be given below. 

(1.3.5) Proposition Suppose that (S ,g) is a simple surface of revolution. Then there exists a coordinate 
system (u, 9) on U, with sinu = a(r) and a smooth function f on [—1, 1] such that /(l) = 1, /(— 1) = 1 and 
such that 

g = [f (cosu)] 2 du 2 + sin 2 ud6 2 . 



Proof: First, define 



and 



b(r) 



sin (a(r)) r € [0, r ) 

ir — sin^ 1 (a(r)) re[r a ,L] 



Then define / : (-1, 1) ^ IR by 



M = j (a|[o,r,])- 1 (v / IHg) v€[0,l] 



0] 



Since b{r) = u,c(cosu) — = r we have a(r) = sinu,a'[c(cosu)]dr — cosudu and b'(r)dr — du. 

The smoothness properties of / follow from those of a and from the fact that a(r) and sinu have the same 
qualitative shapes [Besse, loc.cit.]. ■ 

The geometric result is the following [Besse, Theorem 4.11]: 
(1.3.6) Proposition In the above coordinates, the equator 7b(s) = (u(s),9(s)) has the equation 

«( s ) = |, d(s)=s. 

Any other geodesic 7(s) = (u(s), 0(s)) in U is contained between two parallels u — i and u = n — i and the 
angle 9(i) between two consecutive points of contact with these parallels is given by: 



9(i) = sini 
The length of this arc of 7 is given by 

s(i) = I 

J i 



f(cosu) 



sinu(sin 2 u — sin 2 i) 
sin(u)f(cosu) 



—du. 



(sin 2 u — sin 2 i)2 



—du. 
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Figure 3: 



Sketch of Proof: Using the Clairaut integral, the equations of the geodesic have the form: 



( de 

ds 
d0_ 

du 



sin^u 

sin(i) f (cosu) 



sinu(sin 2 u— sin 2 i) 2 
ds sin(u)f(cosu) 

du (sin 2 u- sin 2 i) 2 

The formulae above follow by integration. ■ 
The following is geometrically obvious: 

(1.3.7) Proposition Let Tj C S*S 2 denote an invariant torus with H(I) — 1 and let i(I) be the u- 
coordinate of the extremal parallel closest to N in the projection of Tj to S 2 . Then: 
(i)T E (I)=2s(i(I)); 



(ii) lo e (I) = 2(0(<(I)) - 1); 

(in)r E {I) = lf r ^]{l-^)-Ur 

(iv) u E (I) = -1 + £ 0% a(rr 2 (l - ^)~Ur. 



Proof: 

(i) By definition, lue(I) is the change in angle along the equator between a geodesic 7/ on Tj, starting at 
the equator on a fixed meridian and heading towards the northern hemisphere, and the fixed meridian, upon 
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second intersection with the equator. We can view this arc of the geodesic as consisting of three pieces: 
one from the equator to the northern extremal parallel, one on the 'back-side' between the two extremal 
parallels, and one on the 'front-side' from the southern parallel to the equator. Since the lengths of the 
'front-side' arcs are unchanged by rotation, we can rotate one until the two make up a smooth geodesic arc 
between the parallels. The length of the geodesic arc is therefore twice that of an arc between the parallels, 



(ii) For the angle change: In the same way, the change in 9 along this arc of 7/ is the change in 9 of two arcs 
between the extremal parallels. We subtract 1 since utj measures the addition to one full revolution. 

(iii) Since sin(u) = a(r) and f(cosu)du = dr, we have 



(iv) From Propositions 1.3.2 we have that ^ = j^uje- Since H(Ii,I 2 ) = 1 implies that w\ + F'{I\)uj2 = 0, 
we get from (1.3.1) that lu e = -2irF'{h). ■ 

(1.3.8) Corollary The non- degeneracy condition (1.1.6) is satisfied if V 1E is globally twisted, i.e. if 
lj' m > or lu' m < 0.. 

Proof: In these cases, F is convex (or concave). Since the set {H = 1} is the graph of F in r o , the Gauss 
map (and hence the frequency map) is an embedding. ■ 

Remark Both cases, of concavity and convexity, occur for ellipsoids of revolution, see [Bl]. The separating 
case of the round sphere is of course degenerate. 

1.4 Classical Birkhoff invariants 

The classical Birkhoff normal form of Hamiltonian H near a non-degenerate periodic orbit 7 is a germ of a 
completely integrable system to which H is symplectically equivalent in a 'formal neighborhood' of 7 (see e.g. 
[F.G] for a detailed discussion). In the case at hand, where H is already completely integrable, the Birkhoff 
normal form is simply H itself expressed in terms of action-angle variables. The 'Birkhoff invariants' of the 
Hamiltonian at a torus T] o may then be identified with the germ of H(I\, I2) or of u>i at I = I . 

Since H is homogeneous of degree 1, it is equivalent and somewhat clearer to define the Birkhoff invariants 
after first breaking the JR + symmetry. That is, we would like to introduce a 'base' to the cone r o . The most 
natural one may appear to be the energy level {H = 1}; but for the purpose of calculating wave invariants 
at Tj o it is more convenient to use the tangent line luj o ■ (I — I Q ) = at a point I £ {H = 1}. 

Let us first consider the level set {H — 1} as the transversal. The homogeneous function H is obviously 
determined by the curve H(Ii,I 2 ) = 1 whose equation is given by (1.3.1). Hence we can define the Birkhoff 
invariants at a torus Tj o to be the Taylor coefficients of the function F' at I Q £ {H = 1}. By (1.3.7 (iv)) it 
is equivalent to put: 

(1.4.1) Definition The (first) Birkhoff invariants of H at an invariant torus T L with 1 £ {H = 1} are the 
Taylor coefficients of ue at 1 in the coordinate I\. 



i.e. t e (I) = 2s{i(I)). 
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Secondly, let us consider the tangent lines as transversals: We fix a point 1° E {H = 1}, let ui° denote 
the common frequency vector of the ray of torii IR + T/° and put 

2 

(1.4.2) :=5>£/ fe . 

fe=i 

The equation of the tangent line to {H = 1} at 1° in the action cone T is then given by lo° ■ I = 1. Note 
that I ■ uj° \s homogeneous of degree 1 and hence equals H(I) along the ray H/°; consequently it is elliptic 
(non- vanishing) in a conic neighborhood W C r o of it. The conic neighborhood will be parametrized in the 
following way: we fix a basis (i.e. a non-zero vector) v of the line I ■ uj° = 0, and define the map 

(1.4.3a) M->p(I° + tv), £€(-e,e). 

For sufficiently small e, this map sweeps out a conic neighborhood W a of 1° with inverse given by 
(1.4.36) p = w°-I, i>, : /;• 

Since H{h,I 2 ) = (lo° ■ I) H(^j, £j) and since ^) e {w° ■ I = 1} we may write 

(1.4.4a) H{h,I 2 )=ph I o(£) 
where hio is the function on W n {w° -1—1} defined by 
(1.4.46) /i JO (0 :=H(r + Zv). 

The C°° Taylor expansion of around ^ = is then a symplectic invariant of H. 

(1.4.5) Definition The second (tangential) classical Birkhoff invariants of H associated to the periodic 
torus Tio are the Taylor coefficients hfo (0) . 

In the real analytic case the Taylor coefficients determine hjo and hence H by homogeneity. It is more 
or less obvious that the first and second Birkhoff invariants also carry the same information in the C°° case. 
To be sure, we prove: 

(1.4.6) Proposition The first Birkhoff invariants canonically determine the tangential Birkhoff invariants 
and vice versa. 

Proof: 

By definition 

h IO (0=H(I°+^v). 

Therefore, 

H(h I o(Z)-\l° + Zv)) = l 

or equivalently 

(o - 1 ( j 2 ° + & 2 ) = F(h IO (o - 1 (^i + n ) ) . 
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Writing u — , this says 

F(u) = I 2 ° + jv 2 
u If + &i ' 

Hence the knowledge of Taylor coefficients of F is equivalent to knowledge of the Taylor coefficients of hjo . 




Figure 4: 

We may reformulate the non-degeneracy condition (0.1 (iii)) in terms of where Tj° is the meridian 
torus: 

(1.4.7) Proposition (S 2 ,g) satisfies the non-degeneracy condition (0.1 (Hi)) as long as a := hj o (0) ^ 0. 

Proof: By definition, = H{I° + £v). At the meridian torus, v = (1,0) so = H(I° + (g,Q)) and 

Ko(0 = wiHil + K,0)) = wi(I° + K,0). Hence 

^(0) = ^-^(I°). 

Since wi(/°) = it follows that 

Hence, the condition h"„(0) ^ is equivalent to the condition that ■^■lum(I ) =/= 0. The equivalence of this 
to (0.1 (iii)) is proved in Propositions (1.3.4) and again in Corollary (1.3.8). ■ 
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1.5 Symplectic conjugacy of geodesic flows 

The Birkhoff normal form of a Hamiltonian H at a closed orbit (or family of closed orbits) is a symplectic 
conjugacy invariant of H in a neighborhood of the orbit (s). Hence the global Birkhoff normal form H(Ii, I2) 
of a completely integrable Hamiltonian is a symplectic conjugacy invariant. The purpose of this section is 
to show that it is a complete conjugacy invariant. We begin by showing that the homogeneous Hamiltonian 
torus actions commuting with geodesic flows of simple surfaces of revolution are all symplectically equivalent: 

(1.5.1) Proposition Suppose (S 2 ,gi) and (S* 2 ,^) are smooth surfaces of revolution of simple type and 
let (Ji.,/2) resp. (Ji, J2) denote their global action variables as above. Then there exists a homogeneous 
symplectic diffeomorphism \ '■ T* S 2 — > T* S 2 such that \* Ji = I{. 

Sketch of Proof : Let fa, resp. ipi be the dual angle variables on (S 2 ,gi) resp. (S 12 ,^)- Except on 
the equators of (S 2 ,gi) a point in T*S 2 is uniquely specified by its action-angle coordinates. Define \ to 
be the identity map in action-angle coordinates with respect to the two metrics. It is obvious that x is a 
homogeneous symplectic diffeomorphism on the complement of the equators, so to prove the Corollary it 
suffices to show that \ extends to the equators with this property 

Since \ is homogeneous of degree 1, it is necessary and sufficient to define it on the unit cotangent bundles. 
Moreover, since it commutes with the Hamilton flow of pg on the regular set, its extension must also do so. 
Hence it must be the lift of a diffeomorphism \ on the orbit space O := S*S 2 /S 1 of the rotation action. This 
action is free so the natural projection p : S*S 2 — > O must be diffeomorphic to the standard projection from 
50(3) — > S 2 . The image of the torus foliation defined by level sets of pg is a singular foliation of O formed 
by level sets of the function pg induced by pg on O, and the two singular points are the images of the 
equators 7^. Since pg is, by assumption, a perfect Morse function on O ~ S 2 for each metric, the quotient 
map x on the punctured quotient O — {o ± } of \ extends smoothly to the completion. It follows that \ 
extends smoothly as a rotationally equivariant map on the completion of S*S 2 — {7^} and the homogeneous 
extension must be symplectic. (See [CV.l] for more on the behaviour near the poles). ■ 

Remark: In fact, this proposition can be sharpened to say: there exists only one homogeneous Hamiltonian 
S 1 x S 1 action on T*S 2 — (up to symplectic equivalence). This follows from a homogeneous analogue of 
the Delzant classification of completely integrable torus action on compact Kahler manifolds. We hope to 
report on more general results of this kind at a later time. 

Thus, the torus actions defined above on the cotangent bundles of simple surfaces of revolution are 
always symplectically equivalent. The question arises when the geodesic flows are symplectically equivalent. 
The answer can be given by expressing the norm functions |£| fl of the metrics in terms of the global action 
variables. Before doing so, we note that the action variables are not quite uniquely defined above because 
the choice of generators 7,(6) is not unique. For instance, one might have permuted the roles of N and S. 
Hence we have: 

(1.5.2) Proposition Let (S* 2 ,^) be a simple surfaces of revolution, and let (Ii^h) resp. (Ji, J 2 ) be the 
global action variables defined above. Let = _ffi(7i,/ 2 ) resp. |£| fl2 = ^2(^1, J2) be the expressions of 
the metric norms of gi in terms of action variables. Then: the geodesic flows of(S 2 ,gi) are homogeneously 
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symplectically equivalent if and only if there a linear map A = [ 011 0-12 ) e SL(2, 2Z) such that Hi = 

\ «21 »22 / 

H 2 O A. 

Proof If such a choice exists, then the map \ above obviously defines a symplectic conjugacy. 

Conversely, suppose the geodesic flows are symplectically conjugate by a homogeneous symplectic diffeo- 
morphism x, i-e. x* H 2 (Ji 1 J 2) = Hi (Ii,I 2 ). Then x* Ji are global action variables for the geodesic flow of 
(S 2 , gi). But global action variables are almost unique: they correspond to a trivialization of the lattice bun- 
dle Hi(F b , ZZ) [D][G.S]. Therefore there exists A e SL(2, 7Z) so that X * J = A ■ I. Then Hi(I) = H 2 (A ■ I), 
m 

The following gives a geometric interpretation of the conjugacy condition: 

(1.5.3) Proposition Suppose gi,g 2 € 1Z. Then their geodesic flows are symplectically equivalent if and 
only if their equatorial first return times te and angles lu^ are equal. 

Proof: Suppose first that the flows are symplectically equivalent. After choosing compatible bases for the 
homology we may then assume that the expressions for Hi and H 2 in global action-angle variables are the 
same. Then the frequency maps are the same and hence by Proposition (1.3.3) the equatorial return angles 
are the same. Also, the Poincare maps V 1E are conjugate and hence the equatorial return times are equal. 

Conversely, if the equatorial return times and angles are the same, then the flows have the same frequency 
vectors (as functions of the global action angle variables) and hence have the same global Birkhoff normal 
forms. By Proposition (1.5.2) the flows are symplectically equivalent. ■ 

2 Quantum dynamics and normal form 

It is owing to the following notion that simple surfaces of revolution are so manageable. 

(2.1) Definition The wave group e 4tv ^ of a compact, Riemannian n-manifold (M,g) is quantum torus 
integrable if there exists a unitary Fourier-Integral representation 

r:T n ^ (7(L 2 (M)), f (tl; ..., tn) = e *(*iA+-*»'») 
of the n-torus and a symbol H € 5 fl (K™ — 0) such that \/A = H(Ii, ■ ■ ■ , I n )- 

In the formula above, we follow physics notation in indicating operators (as opposed to their symbols) 
with a 'hat'. Thus, the generators Ij are first order pseudodifferential operators with the property that 
e 2mI i — Cjld for some constant Cj of modulus one. 

Since if is a first order elliptic symbol on 1R™ — it has an asymptotic expansion in homogeneous functions 
of the form: 

(2.2) H ~ Hi + H a + + . . . , Hj (rl) = r^H 3 (I). 

The quantum action operators are uniquely defined up to a choice of basis of H 1 (fii n /2Z n , 2Z), the terms 
Hj are uniquely determined up to the same ambiguity. 
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The principal symbols Ij of the Ij 's generate a classical Hamiltonian torus action, so any quantum torus 
action automatically induces a classical one. Conversely, it is a theorem of Boutet de Monvel-Guillcmin and 
Weinstein that any classical Hamiltonian torus action can be quantized [BM.G, Appendix, Proposition 6.6]. 
Since metrics in 1Z give rise to quantum torus actions, we have, in particular: 

(2.3) Proposition (cf. [CV.l]) Suppose g € 7Z. Then J A g is quantum torus integrable. 

We also observe that the following holds for any Laplacian commuting with a quantum torus action: 

(2.4) Proposition For any A, H D = 0. 



Proof: Since the subprincipal symbol a su b(V A) equals zero, the same is true H := H(Ii,Ii). Now <7 sub (H) 
is invariantly defined (independent of a choice of symplectic coordinates); hence it may be expressed in 
action-angle coordinates in the form 

But the mixed derivative term automatically vanishes since Hi is a function only of the classical action 
variables. ■ 



Let us specialize to the case of wAg with g e 1Z. From the fact that e 3 ' = Cjld for a quantum torus 



A 9 

action, it follows that the joint spectrum of the quantum moment map 

Sp(X) C 2Z 2 n T + 

is the set of integral lattice points, translated by /i, in the closed convex conic subset T C 1R 2 . The vector 
/i = (/ii,yU2) can be identified with the Maslov indices of the homology basis of the invariant torii. In the 
case of V /SJ, n = (0, 1/2) [CV.l]. 

Expressing a/A^ in the form H(Ii,I 2 ) we have that 

Sp{^/A~ g ) = {H(N + n) : N e 2Z 2 n T a }. 
Thus, the eigenvalues of y/~A~l have the form: 

\ N ~H 1 (N + m) + H_!(N + + 

and the wave trace takes the form 



(2.5) Tr ^V^ = < 



The symbol H may be regarded as a global quantum Birkhoff normal form. As in the classical case, it 
has a germ at any periodic torus, so these may be regarded as the microlocal Birkhoff canonical forms. To 
be more precise, we imitate the definition of the classical tangential Birhoff normal forms and write 

(2-6) H j (I 1 ,I 2 ) = (u°.iy Hji-^-^^y^^o.jy hj(0 
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Figure 5: 

where £ is a linear coordinate relative to a vector v generating ui° ■ I = 1. We then Taylor expand hj(£) 
around £ = 0: 

a>0 

(2.7) Definition TVie homogeneous quantum Birkhoff normal form coefficients of H at the periodic torus 
Tjo are the Taylor coefficients h" for j = 1, 0, — 1, . . . . 

(2.8) Proposition Assume that g € 1Z* . Then all of the functions Hj are real analytic in the interior of 
the cone T a and all of the functions hj are real analytic near £ — 0. 

Proof: First consider |£| 9 = Hi(Ii, I2) = H(I\,l2). We know that H is a C°° homogeneous function on r o . 
On the other hand, from the formula I2 — G(|£| g ,/i) we see that I2 is a real analytic function of |£ 3 |,ii. 
Since G is the inverse function to H with respect to the first variable, H must also be a real analytic function 
of/1,/2. 

Next recall that I2 is defined (up to a smoothing term) as the function G(V~K,Ii) which has principal 
symbol G(|£| s ,Ii) and which has integral spectrum. More precisely, one begins with G(v / A, h), which has 
the property that exp(2-KiG{\f~K, I\) + fa) = I + K with K of order -1. One defines R = ^pLog(I + K) and 
puts <5(\/A, h) = G( v / A, h) + fa + Rj- 

Since G is a real analytic function and we only apply the holomorphic functional calculus in the steps of 
the construction, it follows that G is an analytic function of (\/A, 
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Since 

h = G(\/A,7i) has the inverse function \/~A = H(I\,l2), it follows again by the inverse function 
theorem for analytic functions that H is a real analytic function. 
The real analyticity of the hj's follows from that of the Hj's. m 

3 Wave invariants as non- commutative residues 

To relate the wave invariants to the coefficient of the normal form, it will be helpful (as in [G.l] [Z.1,2]) to 
express the wave invariants as non-commuative residues of the wave operator and its time-derivatives. Let 
us recollect how this goes. 

The non-commutative residue of a Fourier Integral operator is an extension of the well-known non- 
commutative residue of a pseudodiffcrential operator A on a compact manifold M , defined by 

rcs(A) = 2 Res s=0 Q{s, A) 

where 

C(s,A) = Tr AA- S / 2 (Re s » 0) . 

Here, A is a Laplacian (or any positive elliptic operator) on M . From work of Seeley, Wodzicki and Guillcmin, 
one knows that ((s, A) is holomorphic in Re s > \ dim M+ ord(A) and admits a meromorphic continuation 
to C, with simple poles at s = dimM+ord (A) — k (k = 0, 1, 2, . . .). The residue at s = has a number of 
remarkable properties (not shared by the residues of the other poles): 

- ies(AB) = res(BA), i.e. res is a trace on the algebra \I/*(M) of pseudodiffcrential operators over M; 

- res (A) is independent of the choice of \/A; 

- there is a local formula for the residue, res(^4) = (27r)~™ J S , M a-n(x, £)indx A d£, 

— oo 

where a_„ is the term of degree (— n) in the complete symbol expansion a <~ ^ aj for A; dx A g?£ is the 

m 

canonical symplcctic volume measure on T*M. 

These results may be extended to Fourier integral operators as follows: Let A be a Fourier Integral 
operator in I m (M x M, A) for some homogeneous canonical relation A C T*(M x M)\0 and m € 7Z. Below, 
diag(X x X) denotes the diagonal in X x X. Below we will sketch a proof of the following: 

(3.1) Theorem Suppose A and diag(T* M x T* M) intersect cleanly. Then ((s, A) = Tr AA~ S / 2 (Res » 
0) has a meromorphic continuation to C, with simple poles at s = m + 1 + — j, where eo = dim An 
diag{S*M x S*M), and j = 0, 1, 2, ... . 

The clean intersection hypothesis above is that An diag(T*M x T*M)\0 is a clean intersection. It is 
satisfied in the case where A = C t if and only if the fixed point set of G" is clean. Hence, Theorem (3.1) 
implies: 

(3.2) Corollary Let ((s,t) = Tr U(t)A- s ^ 2 . 
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If the fixed point set of G 1 is clean, then ((-,t) has a meromorphic continuation to C, with simple poles 

at^i+ ^y-' -j o =0,1,..., 2). ' 

Here, S Fix(G*) is the set of unit vectors in the fixed point set of G*. In the case of a completely integrable 
system, S Fix(G L ) is the union of the periodic torii with period L. We assume here, and henceforth, that 
the periodic torii are all clean fixed point sets for G* on S*M. 

The non-commutative residue of the Fourier integral operator A is then defined by: 

(3.3) Definition 

res(A) := Res s=a ((s, A) 



just as in the case of pseudodifferential operators. And just as in that case, res(A) has some remarkable 
properties: 

- it is independent of the choice of A; 

- if cither A or B is associated to a local canonical graph, then ves(AB) = tcs(BA) 

- there is a local formula for ves(A). 

The basic properties of res(A) may be deduced from a singularity analysis of the closely related distribu- 
tion trace S(t, A) := TrAU(t) (cf. [Z.4]). Under the cleanliness assumption above, S(t,A) is a Lagrangean 
distribution on 1R with singularities at the set of 'sojourn times', 

Sr = {T: 3(x,£) G S*M : (x, £, G T (x, £)) e A}. 

For a given sojourn time T, the corresponding set Wt — {(x,0 '■ (x,^,G t (x,^)) <G A} of sojourn rays fills 
out a submanifold of S* M of some dimension &r- With N(X, A) =J2 /\:<\(A<Pj,<Pj), we have 



(3.4) Proposition [Z.4, Proposition 1.10]. with p e C£° (B) , supp (p) n ST = {0} and 

p = 1 near 0, then 

p * dN(\, A) - G„A m+£ V L a-jX- 1 , 
where C n is a universal constant. The coefficients have the form, 

U>jd\A 

Aa 

where Aa = An diag(S* M x S*M), eo = dim Aa, d\\ is a canonical density on Aa and the functions ujj 
are determined by the j-jet of the (local) complete symbols of AU along Aa- 

Proof of Theorem (3.1). As in the pseudodifferential case [DG, Proposition 2.1], we have 

XrAA - s / 2 = ( XsjCW (-,A)) 

= (x s ,p*dN(-,A)) + entire. 

Hence, 



oo POO 

((s,A) = C n J2<*j X^ +m - s -'d\- 

4— n J 1 



entire 



3=0 
oo 



= C n > t h entire 

"^m+l + ^-( S + j) 
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completing the proof. H 

In particular, if A = U(L), then {0} is a sojourn time if and only if L e Lsp(M, g). If L £ Lsp(M, g), 
C(s, L) is regular at 0. 

Now let us return to the case where the geodesic flow is completely integrable. Then the dimension of 
each periodic torus T of period L equals e Q = dimT — n. Hence we have: 

(3.5a) TrU{t) = e a (t) + ^ e T (t) 

T 

where the sum runs over the periodic tori in S*M and where 

(3.56) e r (t) = a T ._n±i (t - L + iOy^ + a T ._n±i +1 (t - L + iQ)~ s ^ L+1 + .... 

More precisely, it takes this form if n is even; if n is odd, the positive powers of (t — L + iO) should be 
multiplied by log(t — L + iO). In the following Corollarly we use the notation — T for the time reverse torus 
crT. 

(3.6) Corollary 

^ n+l | f. 

22 a ±r,-(ii±±)+k = res(VA 2 U(t)\ t=L ) 
± 

Proof. In the notation of (possibly negative) fractional derivatives in t, we have ^U(L) — D^U{t)\ t= L =■ 
The claim follows from the facts that 

D;^ +k (t -L + iO)-^^ = log(t -L + iO) 

and that the non-commutative residue is the coefficient of log(t — L + iO). m 

(3.7) Examples: 

(a) If dimM = 2, then the wave trace expansion at a torus T has the form 

e T (t) = a T _3 (t - L + i0)~i + a T _i(t - L + i0)~i + . . . . 

Hence 

± 

(b) If dimM = 3, the wave trace expansion at T has the form 

e T (t) = a T - 2 (t - L + iOy 2 + a T ,-i(t -L + iO^ 1 + a rfl log{t — L + i0) + .... 

Hence 

]>> ± T,-2 = re S A- 2 e 4Lv ^, a T ,-i - resA^e 1 ^, ]T a±r,o = rese lL ^, . . . . 
± ± 

(3.8) Remark: The wave invariants for a closed geodesic 7 (or periodic torus T L ) are exactly the same 
as for its time reversal, hence the same residue formulae also give the individual wave invariants. For this 
reason it is often not necessary to resolve the ambiguity between a torus and its time reversal. 
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4 Proof of the Main Lemma 

We now prove that the wave trace invariants of a metric g G TZ* determine its quantum normal form H. 
4.1 Wave invariants for simple surfaces of revolution 

We begin by specializing Corollary (3.6) to the case of a surface of revolution (S 2 ,g) in TZ* . Since the length 
spectrum L is simple, there is a unique periodic torus Tl C S*S 2 of each length L € L (up to time reversal). 
By the existence of global action-angle coordinates, it may be expressed in the form Tl = Ti L for a unique 
point II = (Ili,Il2) € V (up to reflection). Let wl denote the frequency vector at II- Then we have 
L ■ lul = Ml where Ml is the vector of winding numbers of the periodic orbits on Tl- 

The other periodic torii of period L lie on the rays !R + 7i U IR + (t(7l). Their action coordinates lie 
on the rays M + Il U 1R + <j(Il) and they have the same frequency vector, u>l, as for II- Since the wave 
invariants at Tl depend only on the microlocalization of the wave group to a conic neighborhood of TR + Tl , 
we introduce a microfocal cut off operator V>z,(/i, h) with ip homogeneous of degree 0, equal to 1 in a small 
conic neighborhood of the ray M + Il and zero off of a slightly larger conic neighborhood. The singularity of 
TrU(t) at t = L is then the same as the singularity of Tr^pL U{t). 

To calculate the wave trace as a residue, we also introduce the first order pseudodiffcrcntial operator 
u>l ■ I '■= ULih + ^L2-^2- We emphasize that u>l ■ I is a linear combination with constant coefficients of the 
action operators. Since u>i ■ I = H, the principal symbol wl • I of wl • I takes the value 1 at I = II and 
therefore Wl • I is elliptic in a conic neighborhood of 71 . We will use it as the gauging elliptic operator in 
the residue formula for the wave invariants. 

(4.1.1) Proposition Let g G TZ* and let Le£. Then we have: 

E a ±T L ,-i+k = Res s =o I MI + ^ M ^e-^ {I+ ^{H{I + ^ +k {^L ■ (I + n))-'dl 
± "' r 

where as above Ml is the vector of winding numbers of Tl . 

Proof: Since \fK = 77(/i,/ 2 ) and since ipL is a function of the action operators, we have by Corollary (3.6) 
that 

£ a ±T„-i+* = i?es s=° J2 + ti)e^ N+ ^(H(N + M ))^ +fe K • (N + M ))" s . 

± Nezz 2 

We then apply the Poisson summation formula for Re s » to replace the sum over N e 2Z 2 by 

Res s =o ^2 Jl,mM s ) 

where 

Jlm.M ■= MI + li)e-^ M ^e iLk ^\H{I + ^))^ +fe K ■ (J + ^))- s dI. 

It follows from Theorem 3.1 and by simplicity of the length spectrum that only the term with M — Ml 
has a pole at s = 0. This can be seen more directly from the fact that only in this term does the phase 
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— (M,I) + LH(T) have a critical point, since M = LViH(Im) implies that the torus with actions Im is 
periodic of period L (cf. §1.2). ■ 

To calculate the residue of the integral Jz,fc(s) '■= JL,M L .ki s ), we rewrite the integrals in terms of the 
(p, £) coordinates introduced in (1.4.3) in a conic neighborhood of the point II € {H = 1}: Recall that we 
parametrize the tangent line T Il {{H = 1}) by 

where v is a non-zero vector along the line Ml ■ I — and parametrize a conic neighborhood of M + Il by 

(p,0 eM+ xM^ P (I L +£v). 

Since I L is fixed, we abbreviate hj L by h in the next proposition. We also denote the cutoff in these 
coordinates by V>l(0- 

(4.1.2) Proposition With the above notation, a Th _i +k equals the term of order ( o~ 1 ~ fc in the asymptotic 
expansion of the integral 

J L k(p) := L - 2+s - k+ h l < M ^ p? f ^ L (C)e VQ ^e W3 e l Sr. 1 L ' c+1 ^' c ' l -(«)[/ l (e) + f;^+V- £ - 1 ^(0]"^ + ' c rfe 

Here, a = h"(0), and gs is the third order remainder in the Taylor expansion of h(£) at £ = 0. 

Proof We first change variables I + p — * I in the expression for JL,M L ,k(s) in the preceding Lemma and 
then further change variables to (p := % • J, £). Thus we get 

J Lk ( s ) = e i < M '"> [°° f M0e~ lp{M ' {lL+ ^e lLlph(i)+p " h - l{0+ ^ 

Jo JTR. 

Taylor expanding = /i(0) + h'(0)£_ + \h" '(0)£ 2 + g 3 (£) and using that h(0) = H(I L ) = 1, that h'(0) = 
uj Il ■ v = and that (M, I L + £v) = (M, I L ) = L by (1.2.1a-b), we get 

j Lfe ( s ) = e *<^> P f ^(Oe^'w^V^ 3 ^^ 

Jo JjR. 

Now change variables again, p — > L -1 /? and pull the factor pi -1 in front of /i(£) in the bracketed 
expression outside the d£ -integral. We get: 

/*°° /■ 

J Lfe ( s ) = e «M,,.> £ -2+.-fc+i / { / ^ L ^y P h''me e i P m{O e iL[L P -'h.^)+...] 

Jo JjR. 

[h(0 + p- 2 L 2 h^(0 + . . .}- i+k d^}p- s+1+k -^dp. 

The pole at s = is produced by the terms of order p~ 1 in the d£ — dp- integrals, hence by the terms of order 
— k — 1 — | in the asymptotic expansion of the d£-integral. ■ 
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To determine the terms of order p k 1 in JLk(p) we apply the method of stationary phase. Cancelling 
factors of p? we get: 

(4.1.3) Corollary In the above notation, a TL _i +k equals the term of order p~ k ~ x in the asymptotic 
expansion of 

oo oo 
_J =L -2-k+i e r{M^ ^ (2z p)-™ a -™a| m { e ^(« e l ^i iJ + 1 ^' l -(«[/ l (e)+^^+V- £ - 1 / l -£(e)]-^ + ' £ }|c = 



4.2 The meridian torus 

We now want to analyse the terms of the expansion coming from the ray of meridian torii with I\ = 0. 
The special property of the meridian torii is that they are invariant under the canonical involution a(x, £) = 
(x, — £) of T*S 2 — 0. As the following proposition shows, this leads to a useful symmetry property of the 
quantum normal form. As above, the notation h(£), etc. refer to the functions hj L (etc.) where II is 

the point on {H = 1} corresponding to the meridian torus in S*S 2 . 

(4.2.1) Proposition The complete symbol of H is invariant under a. Hence hj is even for all j. 

Proof Let C denote the operator of complex conjugation: Cip — 4>. Since \/A commutes with C, so does 
H(I). Since C is a conjugate-linear involution, this implies that 

H(I) = C^H^C = H(C- 1 IC) 

where the bar denotes complex conjugation. 

We claim that C~ X I\C = —I\ and that C~ x IiC = I\. Moreover, that I\ o o = o a = I 2 - The 

statements regarding I\ , I\ arc obvious since I\ — -J^q chan ges sign under complex conjugation and since 
o~ c -if iC = I\ ■ o~. This latter also follows from the fact 7i(x,£) — X2£,\ ~ xi^2- 

Regarding I2 we note that its cr-invariance follows immediately from the explicit formula (1.3.1). The 
non-obvious claim is that I2 is invariant under C. But from the invariance of the principal symbol we have 

C- 1 / 2 C = / 2 +if 1 

where K\ is of order 0. Since I 2 is a function of (7i,\/A), it is clear that C must take joint (h,l2) 
eigenfunctions into joint eigenfunctions. Hence C determines an involution on the joint spectrum {(n, k+ i) : 
\n\ < 2k + 1, k > 0} which, we recall, is simple. Consequently, the involution (still denoted C) must take the 
form 

C : (n, k + -) -» (-n, k + - + f(n, k)) 

where / is a bounded function. Moreover, since f(n,k) — (l2<fi n ,k, <t>n,k) and since the § n ,k are quasi-modes 
associated to Bohr-Sommerfeld-Maslov torii [CV.3], it is clear that /(n, k + i) must be a polyhomogeneous 
function of order on r o . Since it is also integral- valued on the semi-lattice of joint spectral points, it must 
be constant. Additionally it must satisfy the involution condition, and one sees that the constant must be 0. 
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Returning to H{h,I 2 ) we have H{h,I 2 ) = C- 1 H{I l ,I 2 )C = H{C- x hC, C~ X I 2 C) = H{-I u h) so that 
H(a, b) = H(—a, b). We then observe that H is a real function (at least modulo terms of order — oo.) To see 
this, we recall that the eigenvalues H(n,k + 5) are real and argue by induction on the symbol expansion 
using that 

H{n, k+^)= H(n, k + ^) + if_i(n, k + ^) + . . . e H. 

First, the principal symbol H is real so we may drop it from the expansion without affecting reality. Assuming 
that H = Hi, H—i, . . . , H-k are real, we may drop them all and get that the tail sum is real. Since it is 
dominated by H-/.-1, this function must be real for all points (n, k + 5) sufficiently far from 0. But by 
homogeneity, H_k-i is then real on the projection of these lattice points to H = 1. The projected points 
form a dense set by Proposition 1.2.4 and hence H_k-i is everywhere real. We conclude that H is cr-invariant. 

Now consider the ray of meridian torii, or more precisely the ray !R + (0, 1) in the action cone T. We note 
that this ray is invariant under the involution or(a, b) := (—a,b) of T, and moreover the meridian torus is 
invariant under a (it 'rotates' the torus by angle n.) From the above, the level set {H = 1} is invariant under 
err an d hence the tangent line at (0,1) is invariant. Evidently it is horizontal in the (7i, I 2 )-plane and err 
restricts to it to the map £ — ► — £. Since the complete symbol of H is or -invariant, the h^s must be even. ■ 

We now go back to the wave invariants associated to the meridian torus and its iterates. Let us write 
the amplitude for the kth wave invariant of the iterate of length L, namely 

(4.2.2) ^((.pl^e'^^'^TO+fl^Mfl]-^ 

1=1 

in the form 

A Lk {i, p) := A Lko (0 + P _1 ^Lfei(0 + • • • • 
Thus, a TL _i +k is the term of order p _fe_1 in 

j —k—1 00 00 

(4.2.3) -=== V V a- m (2ip)- m - e dl m [e i ^A Lk( ]. 

V27riLa m=0i=0 

Expanding the derivatives and using that <?3(£) is even and of order 0(£ 4 ), we may rewrite (4.2.3) in the 
form 

j _ k _i 00 00 m 

(4.2.4) J_]r££ £ E {C mli a-™p-™-^(d?ig 3 ...d?« 93 )d? m -»A LM }. 

Here, the multi-index i-sum runs over q-tuplets with i n > 2, |i| = ^i n = j < m with q < |m, 2q < \i\ = j < 
m. In order that —to — I + q = — k — 1 it is necessary that 

TO<2(fc + l), £<k + l. 

These restrictions follow from the fact that it takes at least 4 derivatives of gz to make a non-zero contribution. 
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4.3 A collection of formulae 



For future reference we assemble some notation and identities regarding the coefficients A Lke and their 
relations to the wave invariants. 

(4.3.1) Notation 



(1) Define: F n (h_ u . . . , h- n ) by [h(£) + ET-i L i+1 p~ 

MO" 4 [1 + £n=2 L n p- n F n {h~\h^ . . . , h- n )] 




(2) Define: G n (L,h-i, . . . ,h-„) by e '^=i p 3 


= 1 + E~=i L n p- n G n (L, h- u fc_„). 


(3) Define F nmi by F„(/j-\ . . . , h- n ) = £„=i /i" 


m Ei=(j 1 ,..,j m ):|i|=n Fnmjh-j, . . . h-j m ] 


(4) Define G nmj by G„(L, . . . , h- n ) = Em=i ^ m 


iT,j=(j u ...,j m y.\j\=n Gnmjh-j, ■ ■ ■ /l-jj 



(4.3.2) Identities 

(1) A LOn = L n h~i[F n + G n + E l+3 =n FiGj] 

(2) A Lk = [h(£) + TZi L^p-^h-t(0] k A LQ . 

(3) A LOn = i"/l" 3 ELl ^2j=(j 1 ,...,j m ):\j\=ni F rimjh~ m + G nm jL m }h^ n . . . h-j m 

+ J2a+b=n Emi=l Em 2 =l Si=(ii,...,» mi ):|»|=o '52j=(ji,...,jm 2 ):\j\=b Fam 1 iGbm 2 jh-i 1 ■ ■ ■ ft-i mi ^-ji • • • ^-jm 2 



(4) = feA L(fc _ 1)£ + Ei 7 -=i+,-+i=<, 7 ->i ^ + 1 fe-jA £ ~ 



The coefficients F nm j,G nm j are universal and hence up to the prefactor of .A^jy is a (non- 

homogeneous) polynomial of degree j in L, in ft -1 , and in the /i_/s with universal coefficients. In particular, 
the first few ALfc/s are given by: 

A Lko = h(£)-? +k , A Lkl = fc(0-' +fc L 2 fc-i, A ife2 = MO - *"^ 3 ^ + i 4 ^! + (fc - ^ 2 ^]- 
It follows that the principal wave invariant is given by 

a TL _a := CL = - 1 e'< M ^> 
' 2 V2iriaL 

and that the higher wave invariants a Th _i +k are given by cl times polynomials in L and in the derivatives 
of h, ft -1 , . . . at £ = 0. For instance, the subprincipal wave invariant in dimension 2 is given given in 
terms of universal coefficients C^ k by: 

a T L ,-i — c L[C ooi d^g 3 A Loa + C' ow A LO i + C 002 d^A L00 ] 
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which is easily seen to equal 

c L [C OM dlh(Ok=o + Coi £ 2 ^-i(0) + Coo2^MO"*l*=o]- 

Above the indices in the coefficients Ck m j 1 j 2 ... have the following meaning: k corresponds of course to 
the k index; m is the power of L and j n are the jet orders of h- n , with with exception of j a which is the jet 
order of h. 

(4.3.3) Proposition We have 

a L ^i +k = c L P k (L, h^(0), h^ 2k+i H0), h-!(p), h { ™\0), h- k {0), h 2 _ k (0), ft_ fc _i(0)) 
where P k is a polynomial with the following properties: 

(i) It involves only the first 2k + 4 Taylor coefficients of h at 0, the first 2k of h-\,..., the first 2k + 2 — 2n 
of h- n the first 2 of h- k and the Oth of h—k—i- 

(ii) It is of degree 1 in the variables h- k -i(0), h^ k (0), . . . , ft/ 2fc + 4 ) (0), and each occurs in precisely one term. 
(Hi) The L- order of the monomials containing these terms is respectively L k+2 , L k+1 , . . . ,L°. 

Proof: To prove these claims we combine the formula for a TL _ k+ i in (4.2.4) with the formulae for the Alh 
given in (4.3.2). Since we may replace the factors of g 3 by h in (4.2.4) and since (4.3.2) expresses A Lke as 
h~i times a polynomial in L and in the /i_j's, it is clear that o Tl is given by a polynomial in the data 

stated above. It remains to prove that the polynomial has the properties claimed in (i)-(iii). 
We prove these claims by proving the stronger statement that 

= CL ® k ( L > MO -i , MO, fc (2) (0). • • • - ^ (2fc+4) (0), fc-i(0). • ■ • - fc?i fc) (0), ■ ■ ■ . fc-fc(O), ft-fc(O), ^-fc-i(0))l£=o 

where Q k has the properties (i)-(iii) for variable £. 

The proof is by induction on k. The properties are visibly true for the principal and subprincipal wave 
invariants. Assume then that they are correct for k < N and consider how things change as N — > N + 1. 
First, the amplitude A L r N+1 \ is given by [h + L 2 p~ 2 h_i + . . ]A LN . Second, we are looking at the term of 
order p~ N ~ 2 rather than p-"- 1 in the asymptotic series (4.2.4). 

In going one further order into the asymptotic series, only two new things happen: 

• The term h-N-2 appears for the first time, arising from the linear term in the Taylor expansion of the 
exponential in (4.2.2). The linear term in the binomial expansion of the power in (4.2.2) does not contribute 
at this stage because its p _1 factor has one higher power. 

• Two additional derivatives are allowed to fall on the previous hj's. A priori, there could be between two 
and six additional derivatives in a method of stationary phase expansion. However, the cases of three to six 
derivatives do not contribute any new data. Indeed, the cases of three-four derivatives occur when m goes 
up by one. One has to remove the extra factor of p^ 1 by applying at least one derivative to the e VS3 factor. 
But in fact all four have to be applied to g% to make a non-zero contribution, and hence no derivatives are 
left to apply to the h-e's. In the case of five or six derivatives, where m goes up to two, one needs to remove 
two extra factors of p^ 1 by applying derivatives to the e V93 factor, and there is no non-zero way to do this. 
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Claim (i) follows immediately from these observations. 

To prove (ii) we note that by the induction hypothesis, /i_jv-i(£)i ■ • ■ j h( 2N+4 \£) occur linearly in Pjy(£). 
Since it requires both of the two new derivatives to fall on these factors to produce h^ N _ 1 etc. in P/v+i(£), 
these factors will also occur to order 1. As for h-N-2, we observed above that it comes only from the linear 
term in the exponential in (4.2.2) and hence it too appears to order 1. 

The proof of (iii) follows the proof of (ii). Since the terms in the exponent of the exponential factor 
in (4.2.2) have the form lJ +1 p~i , the new h-N-2(0 term has the coefficient L N+3 . Similarly, the other 
terms under discussion, e.g. h^^_ 2+r (^), originated as h-N-2+r at the iV + 1 — rth stage with the factor 
of L N+3 ~ r . We would like to show that they remain with just this factor of L as we move inductively up to 
the N + 1st stage. To this end, we note that in order to produce the term h^_ 2 _ r at the (N+l)st stage, 
2r derivatives must fall on the original h-N-2+r- But there are exactly r stages intermediate between the 
N + 1 — rth stage and the N + 1st stage and at each stage at most two new derivatives can fall on a factor. 
Hence, each of the two new derivatives at each stage must fall on the factor of concern. 

Let us also consider what can be multiplied against the original h-N-2+r in the course of produc- 
ing ^7v-2+r(£)- We- observe that each new pair of derivatives is accompanied by a factor of p^ 1 . Since 
/i_Ar_2+r(£), at the N + 1 — rth stage, comes with the factor p~ N ~ 2+r , the 2r further derivatives will bring 
its p _1 -ordcr up to p~ N ~ 2 . Hence no other factors of p^ 1 could have fallen on this factor. Therefore in the 
repeated multiplications by [h + L 2 p~ 2 h-\ + ...], only the repeated choice of the term h, of order in p -1 , 
can have given rise to the term h^_ 2+r . It follows that it retains its original L-order, namely L N+3 ~ r . m 

4.4 Completion of Proof of Main Lemma 

Our purpose is now to show that by using the joint p, L — > oo asymptotics one can recover the complete 
Taylor expansions of all the hj's from the wave invariants of the meridian torus and its iterates. That is, to 
complete the 

Proof of the Main Lemma: We will prove by induction that from the wave invariants dT pL -±+k with 
k < N and all p <= N, we can determine the 2k + 4 -jet of h, the 2fc-jet of h-i,..., the 2k + 2 — 2n-jet of h- n 
at £ = 0. (for n < k + 1.) 

We note that, unlike in the non-degenerate case, the principal wave invariant determines the 2-jet of h, 
and gives no new information under iteration. We therefore begin the induction with the subprincipal wave 
invariant a T _i . 

From the explicit formula above for the subprincipal wave invariant a L _i it is evident that the 2-jet 
term in h is old information, while the other two terms differ in the power of L. Hence they decouple under 
iteration L — > pL and we can determine the 4-jet of h and the 0-jet of h-i from the first two wave invariants. 
Hence the induction hypothesis is correct at the first stage. 

Assume the induction hypothesis for the k— 1st stage. Then the only new information at the kth stage is 
that contained in the terms h-k-i(0), ft^(0), . . . , /i( 2fc + 4 ) (0). By the proposition above, they occur linearly 
in the monomials containing them and the monomials have also the factor of L to the powers L k+1 , L k , . . . . 
Hence the terms containing the new data decouple as p — > oo, and the new data can be determined as stated. 
This completes the inductive argument. 

It follows that we can determine the full Taylor expansions of the hj's at £ = 0. Since they are real 
analytic they are completely determined. Then from the homogeneity of Hj, we can determine Hj from hj 
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and hence the entire function H is determined. ■ 



Remark In the above argument, we are able to drop the terms involving only previously known data 
because of the universal nature of the polynomials Pf.. This universality depends on the fact that we are 
only comparing wave data for quantum torus integrable Laplacians. 

5 Proofs of Corollaries 1 and 2 

We now show that the symplectic equivalence class of a metric in TV is spectrally determined among metrics 



(5.1) Proof of Corollary 1 By the Main Lemma, the principal symbol H g (Ii,l2) as a function of the 
action variables is spectrally determined for a metric g in TV. By Propositions (1.5.2-3), H g determines the 
geodesic flow up to symplectic equivalence among metrics in TV . m 

Remark The Corollary could also be proved by noting that the first return time te is spectrally determined. 
But geodesic flows of in the class TV are symplectically equivalent if and only if they have the same first 
return times te- (See also [C.K] for an equivalent statement.) Also, it should be noted that the first return 
time could be determined from the wave invariants at the equator; hence Corollary 1 would also follow from 
Guillcmin's inverse result for non-degenerate elliptic closed geodesies. 

(5.2) Proof of Corollary 2: Since te(I) is spectrally determined, so is the function s(i(I)) of Proposition 
1.3.6. It is given by 

s(i(I)) = / f(cosu)sin(u)(sin 2 u~sin 2 (i(I))) + 2 du 
Ji{i) 

or, putting x — cosu for u € [0, -|] and — x = cosu for u G [§, t], by 



determines [f{x) + f(—x)]. But [f(x) + f(—x)] is smooth and even so may be written as g(x 2 ) for a smooth 
g; changing variables y — x 2 , and v — u 2 we get 



Thus T is a standard Abel transform, which is well-known to be invertible. It follows that g(y)y 2 is 
spectrally determined; hence so is the even part of /. ■ 



in this class. 




It therefore suffices to show that 
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6 Proof of Final Lemma 



To complete the proof of the Theorem, we need to show that H determines a metric in 1Z. It is plausible 
that this can be done for the following reason: The spectrum of A is the set of values {H(n, k + \) 2 } of H 2 
on the integer points of the action cone V. On the other hand since I\ = the set of values H(n, k + |) 
for fixed n, is just the spectrum {E n k} of the the singular Sturm-Liouville operator 

d 2 , _ _ ^ , n 2 _ 2a(r)a"{r) - (a'(r)) 2 



L n = -{—) +q n (r), q n {r) := q{r) + q(r) = - 



dr ' ' a(r) 2 a(r) 2 

obtained by separating variables in A, fixing the eigenvalue of ^ equal to n, and putting the radial part 
in normal form. Hence, from the coefficients of H we can determine Spec(L n ) for each n. That is, from 
Spec(A) we have determined the joint spectrum of (A, J^). It would thus remain to show that the metric 
can be determined from this joint spectrum, a much more elementary inverse result which has been stated 
several times in the literature ([Kac],[B] [Gur]). Since these prior discussions seem to us somewhat sketchy and 
incomplete, we give a self-contained proof below which was found before we were aware of these references. 



Proof of Final Lemma 



The proof is basically to write down explicit expressions for H and H-\ in terms of the metric (i.e. in 
terms of a(r)) and then to invert the expressions to determine a{r). The first step is therefore to construct 
an initial part of the quantum normal form explicitly from the metric. Up to now, we only know that a 
polyhomogencous normal form exists. 

To be sure, the principal term H has already been implicitly expressed in terms of the metric: Knowledge 
of H is equivalent to knowledge of the level set {H =1} and hence to knowledge of the function F(Ii) in 
(1.3.1). Unfortunately, we have seen in Corollary (5.2) that H only determines the 'even part' of the metric. 
Hence we need to determine at least one of the subprincipal terms H-j. It turns out that only is needed 
in addition to H to determine g. Since the calculation of £/_i requires a new calculation of H, we start 
calculating both from scratch. 

To determine H and H-i in terms of g, we are going to study the spectral asymptotics of VA = H(I\,l2) 
along 'rays of representations' of the quantum torus action, i.e. along multiples of a given lattice point 
(n a , k ). Such rays are the quantum analogue of rays IR + T/ C T* S 2 — thru invariant tori and are basic to 
homogeneous quantization theory [G.S.3]. The basic idea is that the lattice points (n Q , k + parametrize 
tori T no ,k satisfying the Bohr-Sommerfeld quantization condition. To each such quantizable torus one can 
construct a joint cigenfunction <f>n ,k °f (I11I2) by the WKB method. The {<j> n .k} are eigenfunctions of A 
with complete asymptotic expansions along rays. By studying the eigenvalue problem as \(n, k)\ — > 00 wc 
can determine the -ff-j's. 

The WKB method we employ is closely related to the classical WKB method for constructing quasi- 
modes (cf. [CV.3] and the Appendix), except that we have an internal rather than an external Planck 
constant. Let us recall the relevant terminology and notation from the latter case. For each torus TV, we 
denote by 0^(Tj,A) the space of oscillatory integrals (semi-classical Lagrangean distributions) associated 
to Tj, with semi- classical parameters {k m } in a set Ac H + to be specified by the quantization condition. 
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Such oscillatory integrals have the form 

(6.1) u(x,k) = k% + ^J2 I e ikMx >Va{x,£,k)dZ 

where the projection of Tj is covered by open sets lie, where T*(Ue) ~ W{ x V{ and where the phase <t>i{x, £), 
with (x,£) G lie x Ve, parametrizes a part of Tj. The amplitude is a classical symbol in k of order \x. For 
further details, we refer to [C.V.3, §8]. 

As discussed above, the torii Tj project 2 — 1 to the annuli r + (I) < r < r_(7) in S 2 and have fold 
singularities along the extremal parallels. Away from the parallels, an associated quasi-mode is given by 
a sum of two simple WKB functions a±(r, 6){x)e lks ^ ± ( r,e \ Since the actual A- eigenfunctions <f> no ,k are 
quasi-modes attached to the Bohr-Sommerfeld torii T n<j k , they have such a form modulo |(n , fc )|~°°. And, 
since 4> n „, k„ is an exact J^-eigenfunction, its phases must take the form (j>e(r,0) = n o + S no k (r) with 
^-independent amplitudes in polar coordinates. It follows that <t>n ,k (f,^) has the form e m ° f( no ,k ){ r )i 
where f( no ,k )( r ) is an oscillatory integral in the r- variable. It is of course associated to the pushed-forward 
Lagrangean A„ o fc o := p*T, lo k where p : T*S 2 — » T*[0, L] is the projection induced from the map (r, 9) — > r. 
In the case of the meridian torii Tq^ , the pushforward is just a horizontal line p r = C in T*[0, L]. In the 
other cases, the A„ o fc o is a closed curve projecting to an interval [r+(n fc ), r_ (n k )] with fold singularities 
at the turning points (endpoints). The curve is given by an equation of the form H no (r,p r ) = E where 

2 

H n „ (r) = p 2 — -^-pyz is the radial Hamiltonian and where E was the level of the torus. 

The radial part fn ,k (r) of0„ ojfco , is an eigenfunction of the radial operator D no = -(ff-js + ^r 
arising from separating variables in A. Before proceeding, we simplify by conjugating D Uo to the 1/2-density 
radial Laplacian 

~ 1 1 Vb 

D no := a(r)*D no a(r)-* = D 2 r + + W 

a{ry 

where W — a(r)i[— (^) 2 — ^"^] a ( r ) ^- (Note that the volume form dv g — a(r)drd9 of (S 2 ,g) projects 
to a(r)dr.) Thus, we view the radial eigenfunction as having the form g rio ^ (r)\ // adr and apply the WKB 
method to the eigenvalue problem 

(6-2) D no g notko {r) ~ (H(n , k D + i) + #_i(n , k a + i) + . . .) 2 g n ^ ko (r). 

Although the coefficients become singular at r = 0, L, the standard WKB theory applies in the interior 
because it applies to quasi-modes of A on S* 2 . Away from the turning points, the radial part of the 1 /2-density 
eigenfunction therefore has the form 

oo 

(6-3) g no , ko = J2le ±iSn ° k ° {r) J2 a no,k oi m(r)} 

where the phase S Uot k is homogeneous of degree 1 and where the amplitude ot no .k ;m is homogeneous of 
degree — m in (n ,k + 5). These homogeneities replace the powers of k in the non-homogeneous theory 
described above. 

The Bohr-Sommerfeld quantization condition on Tj is that: 

^i=(n ,k ) + ^n , (n,fc)ern^ 2 
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where /j, = (0,2) is the Maslov index (cf. [CV.3, §4]). It is satisfied by T Uo ,k and hence by the radial 
Lagrangean A„ o ,k a ■ The local WKB ansatz (6.3) therefore extends to a quasi- mode of infinite order associated 
to the global A„ oi j. o . Our purpose now is to write down and solve the first two transport equations, which 
are needed to determine H and H-\. For background on the relevant aspects of the WKB method, we refer 
to the Appendix. 

The transport equations are obtained by separating out terms of like order in the asymptotic eigenvalue 
problem 



± m=0 ± m=0 

The leading term, of order 2, is the eikonal equation 

(6.4a) \S' noko {r)\ 2 + H{n , K + l -f 

whose solution is the first order phase function 



(6.46) S noko (r) := J ^ H(n a , k Q f - ^dr. 

The Bohr-Sommerfeld quantization condition on A(„ o fco ) thus reads: 

I no {E):= Area {H Uo < E} = 2n(k + 1) with E = H{n Ql k D + 1). 

The first transport equation, 2a' n k . S' + aS' = 0, is solved by 

1 , i 



(6-5) a (noM {r) = [E noiko - 



a{r) 



2 J 



Hln k +-) 

away from the turning points. Here, E noyko = ° n 2 — — • The solution is of course determined only up to 
a constant, and we have normalized it so that 0L( no ,k o -fi is homogeneous of order 1. 

In the usual way (see the Appendix), we will interpret a not k -o as the coefficient of the ^h„ ~ invariant 
1/2-dcnsity 

vq = [E no ,k r ^]~^Vdr 

on A„ o fe D = {H no (r,p r ) = H(n , k Q + i)}, where ^.H no is the Hamilton vector field of H no . We then re- write 
the higher coefficients a„ ,fc o:m (r) in the form a rio ^ ;m{r)i'o- 

The second transport equation (of order zero) then has the form: 

(6.6) E Hno a notko ,-i = - a ,T„ 1 ,fc„,o(^:) 2Q; «o,fc ,o + W + 2H{n ,k + ^)iJ_i(n , k Q + i), 

where r denotes the local coordinate on h no k obtained by pulling back the base coordinate under the 
projection. The integral of the left hand side over the closed curve A„ o fe o = {H no = H(n ,k + |)} with 
respect to the ^H no -invariant density 

"Lfco.cA = 1 = dr 
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must equal zero. This gives a formula for the first correction to the Bohr-Sommerfeld eigenvalue H(k ,n ): 

I l rr + (n ,k ) d 1 

(6 ' 7) H - l{n °> k ° + 2 )= 2T(n ,k )H(n ,k + l) J rM) ^,0(^)^,0 - ^Wa^dr. 

Here, T(n , k ) denotes the period of the z>h„ flow on A Uo k and r±(n D , k a ) are the turning points. Plugging 

in (6.5) we get, formally, the expression 

(6.8) 

r + (n ,k ) 

dr. 



2T(n , k )H(n , k a + 5 ) J r -(n ,k ) ^E no . ko - ^ 



Actually, the integral in (6.8) diverges at the turning points, and the correct formula is the regularization 
obtained (for instance) by the method of the Maslov canonical operator. For the sake of completeness, we 
provide an exposition of this method in the appendix. Roughly speaking, it regularizes (6.8) by formally 
integrating the 4- derivatives by parts and by moving them outside the integral (in the appropriate way) as 
derivatives in the energy level E. A crucial consequence of the regularization is that only first derivatives of 
a appear in the formulae for H^i(n , k + i). 

Thus the first term of (6.8) is regularized by 



Ci rr + (n ,k ) a ' ( r \2 j ^ 

(6 ' 8 ' lreff) T(n ,k o) H(n , ko + ^J M) " W? T ' 

Here and below d denote (non-zero) constants which we will not need to determine. For the second term 
of (6.8), we note that 

rr+(n ,k ) rr + {n ,k ) rf j ^ rr+(n ,k ) itt'd _i, 

/ Wa nM = ( a a v)— a *dr+ / a a *-—a*dr. 

Jr-(n ,k ) Jr-(n ,k ) « r " r Jr_(n ,fc ) a dr 

After some simplification, this regularizes to: 

C 2 rr+(n ,k ) ffl // V2 1 ^ 



C3 rr+(n ,k ) a /( r )2 1 1 

~T(n ,k )H(n ,k + ±)J r _ {no , ko) W? E ~ W? ] ' r|ls=ls <»»-*»>' 



Now let us return to the inverse problem. We begin from the fact that H(n , 12) is a known function 
of the variable I := I 2 for each n . Its principal symbol H no (I) := H\(n , 1) is then a known function and 
hence the inverse function 

r r+(E) 



Jr-(E) V a ( r ) 
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satisfying H no (I no (E)) = E is a known of E. This of course presupposes that diH no (I) ^ 0, which follows 
from the non-degeneracy assumption (1.1.6). We may write the integral in the form 



(E - x)ldjjL(x) 



where fi is the distribution function of i.e. 



M X) ■= \{r : < x}\ 

with | • | the Lebesgue measure. The above integral is an Abel transform and as mentioned above it is 
invertible. Hence ^ 

dfj,(x) = |— s \~ 1 dx 

is a spectral invariant. Some simplification leads to the conclusion that the function 

J(x):= £ 

r:a(r)=£c 

is known from the spectrum. By the simplicity assumption on a, there are just two solutions of a(r) = x; 
the smaller will be written r_(x) and the larger, r + (x). Thus, the function 

(6-9) J(x) = , * , + 



o'(r_(x))| |o'(r+(x))| 
is a spectral invariant. 

Now let us turn to the H-i expression. Since H(I\, I2) is a spectral invariant, the factors H(n , k a + i) 
and T(n , k ) are spectral invariants. Hence we may remove them from the expression for H—i and still get 
a spectral invariant. For various universal constants C\, C2, C3 it takes the form 

(6.10) [C.dl J L->-(E-^)-*dr + C 2 d E J (E--j)-*dr 
By a change of variables, we may rewrite (6.10) in the form 

(6.11) [dd% J K(x)xi(E-x)~^dx + C 2 d E J K{x)x^{E - x)~^dx 

+C 3 J K(x)x~i(E - x)~ 1 dx}\ E= H( i,,i2) 

where 

(6.12) K(x) = \a'(r-(x))\ + \a'(r + (x))\. 
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All values of E which occur as ratios - U give spectral invariants, so (1.10) (as a function of the variable 
E) is a spectral invariant. 

We now claim that K itself is a spectral invariant. To determine it from (6.11) we rewrite (6.11) in terms 
of the fractional integral operators (cf. [G.Sh, Ch.l §5.5]) 

on the half- line [0, oo]. These operators satisfy 

I a o 1/3 = I a+0 , I_ k = i-j^) k - 
Thus (6.11) equals C{K) where C is the fractional integral operator 

(6.13) C := Ci7_ 3/2 a;5 + C 2 I_ix? + C 3 hx~i. 
To solve for K we apply I_i to CK to get 

(6.14) C[£{x^K{x)) + C' 2 ^K) + C'x^K = I_,CK. 

This equation determines if up to a solution / of the associated homogeneous equation, essentially an Euler 
equation. But also K = on [0, a(r )~ 2 ] Since no homogeneous solution can have this property, K is uniquely 
determined by this boundary condition. 

It follows that both (6.9) and (6.12) are spectral invariants. But from a + b and - + t one can determine 
the pair (a,b). Hence a'(r) is determined from the spectrum. Since a(0) = this determines a and hence 
the surface. ■ 



7 Appendix 

The purpose of this appendix is to give an algorithm for calculating the higher order terms in the quasi- 
classical approximation of eigenvalues for 1 D Schrodingcr operators — — ^3 + V with confining potentials. 
In particular, we carry out the calculation of the h 2 En^ term, which was used in the proof of the Final 
Lemma. 

The algorithm is based on the Maslov method of canonical operators. Expositions and refinements of this 
method can be found, among other places, in Maslov's book [M], in the article of Colin de Verdiere [CV.3] 
(and in its references), and in the recent book of Bates- Weinstein [B.W]. Although these references explain 
the construction of the canonical operator and prove the existence of complete quasi-classical eigenvalue 
expansions, they do not generally go on to describe the calculation of the terms. An exception is the original 
book of Maslov [M], which does calculate the first two or three terms; but the method of canonical operators 
is abandoned at this point in favor of some methods of special functions. As we will show, the canonical 
operator method gives the required corrections quite efficiently. 

The set-up 
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We are concerned with the semi-classical eigenvalue problem: 



\-irfa + V\Mx,h) = E n {h)il> n (x,h) 

(ll>n,lpm) = Smn + 0{h°°) 

E n (h) = Eg\h) + h 2 Ei 2) + h 3 E^ + . . 



The unknown function xp n (x,h) is an oscillatory associated to a Lagrangcan of the form A„ := {H = 
En\h)} where E„\h) is determined by the Bohr-Sommerfeld-Maslov quantization condition: 

Here, /i is the Maslov index of A„; it equals 2 for connected level sets of Hamiltonians of the form H(x, £) = 

e+v(x). 

(k) 

To find the higher order corrections E„ ' we consider the Maslov canonical operator 

U h : S m (A n , Oi ® M) -» O m (TR, A„). 

We follow here the notation and terminology of [BW][CV.3]: S m (A n ,Qi ® M.) is the space of symbolic 
sections of the bundle of 1/2-densities times Maslov factors and m (lR, A„) is the space of oscillatory 
integrals associated to A„. There is a natural symbol map in the reverse direction; any inverse modulo 0(h) 
is a quantization or canonical operator. Its existence is equivalent to the condition that A„ satisfy the BSM 
quantization condition. For background we again refer to [CV.3][B.W]. 

We also denote by Eh the Hamilton vector field of H, by £s H the Lie derivative on any bundle, by s a 
nowhere vanishing section of M., and by p a Cs H -invariant density on A„ (for a fixed n). By surjectivity of 
Uh, we can write an oscillatory integral associated to A„ in the form 

oo 

tp n (x, h)\dx\i - U h [e^ ■ h j fjpi ® s}. 

1=0 

Our aim is to determine the quasi-classical series ^ h? and coefficient functions fj for which the asymp- 
totic eigenvalue problem is solvable. 

We begin by constructing local solutions. Thus we first consider x-projectible pieces of A„ C T*1R: 
pieces which projects regularly from a neighborhood of A € A„ to a neighborhood of x € IR. Restricted 
to densities supported on such pieces, the Maslov canonical operator is truly canonical: if S(x) is a phase 
locally parametrizing A„, then 

oo oo 

U h [e^ ■ Y hi fiP^ ® *] = eiS ■ hia 3 

for some smooth coefficients a,j. We may then substitute this expression into the eigenvalue problem and 
obtain eikonal and transport equations. The eikonal equation (S") 2 + V(x) = E„\h) = has been solved 
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by our choice of phase, so the transport equations become: 

a £zS + 2Va VS = 
ai£?S + 2V &1 ■VS-i£ I a a = 2iE n 2] ] a 
a2^S + 2Va 2 . VS - i^d = 2i[E n 3) a + E n 2) ai ] 

and so on. As is well-known (cf. e.g. [B.W]), these equations may be put into geometric form by observing 
that VS • V is the projection to JR of C~ H and that [aj^S + 2V • S]\dx\ = div(a 2 V S)\dx\ = £ Eh (a 2 \dx\). 
Hence the transport equations become 

£z H (a \dx\i) = 
< C SH ( ai \dx\h) = (2iE n 2) a + i£ s a )\dx\^ 
. £E H (a 2 \dx\-2) = 2i[{Ek 3) a + E n 2) a 1 )+i-£^a 1 )]\dx\i 



The invariant 1/2-density is given by the well-known formula a = (E — V) * , or equivalently by l dx J 2 on 
A„. If we write Oj|dx| 2 = fjpi then / = 1, /j = ^ and the transport equations become 




(2^ 2) +ia" 1 ^a ) 
2i[(^ 3) +£& 2) £)+iaJ 1 -&ai)] 



Here, the expressions in a a , a\, . . . are understood to have been lifted up to A„. 

(k) 

The eigenvalue corrections En are determined by integrating both sides of the transport equations over 
the level {H — En" 1 }. Since the equations are solvable and since the left hand sides will integrate to zero, 
we get 

E " 2) = 2iT(E"») S{ H =E^ a o 1 4s' a o]P 



7?( 3 ) — 1 f rp( 2 )ai , • -1 d 2 „ 1- 



Here, T(E) is the period of 'Eh at level Parametrizing {H = E n ^} as a graph over the x-axis away from 

dx 



the turning points, the invariant density takes the form 7== with £ = E^\ Hence at least formally the 



eigenvalue corrections are given by 

Unfortunately the integral is ill-defined due to the singularities at the turning points. The problem is 
that the Maslov operator cannot be defined near these points as a simple pull-back operator. Rather it 
should be defined as the composition of the Fourier transform with the pull-back operator defined over the 
^-projection. This problem and its solution constitute a key aspect of the Maslov method (in one dimension); 
we refer to [CV.3][B.W] for extended discussions 

The point which is important for us is that the Maslov method gives a regularization of the divergent 
integral. It works in the following way: we introduce a cut-off tps supported away from a ^-neighborhood 
of the turning points x±(E). More precisely we define ipf on A„, equal to one on (25, ^T(E) — 26) resp. 
[\T{E) + 2(5, T(E) - 25) and equal to zero on (T(E) - 5, 5) rcsp. [\T{E) - 5, \T{E) + 5). We then put 

U h (fp ® se^) := htysfp ® *e**) + J h ((l - fa)fp ® se**) 
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where It is the pull-back to IR under the phase parametrization by £ = S'(x) and where Jh is the Fourier 
transform of the ^-parametrization. The notation ipg stands for ipf. For details on Ih, Jh, sec [B.W]. 

Returning to the previous calculation of eigenvalue corrections, we see that what is missing is the cut-off 
tps in the integrals and the contributions from Jh- We wish to avoid confronting the latter. Fortunately, it is 
not necessary to do so: the fact that the eigenvalues are independent of S allows us to determine the Jh (i.e. 
the turning point) contribution indirectly. To see this, we substitute the cut-off into the formula for E n to 
get 

1 r x +( E ) , r p . 

W = ^ l_ (E) (E-V( X ))-^ [{ E-V {x ))-^ X + IIs 

with lis the contribution from Jh . Since the integral is now nicely convergent we can integrate by parts and 
simplify to the form 

1 f*+W V'{xf , A 1 V'(x) , 

The first term tends to 

_1 f_ [ X +W V'(x) 2 

12T(E)dE*J x _ {E) (E-V)^ 5 

as S — > 0. The second expression only involves the Taylor expansion of V near the turning points. Its singular 
part must be cancelled by the singular part of lis, leaving a possible 'residue' as S — > 0. We claim that this 
residue is zero: in fact, this is known to be the case since the first term is well-defined, independent of S, and 
agrees with the formula given in [M]. To give an independent proof that it vanishes, without analysing the 
Jh,-term in detail, we observe that the limit contribution involves only the 2-jet of V at the turning points. 
Hence it must agree with the corresponding expression for a harmonic oscillator at its turning points. But 
no such correction occurs. 
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